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COPYRIGHT (c) 1978, 1980, 1982, 1984 B 
DIGITAL EQUIPMENT CORPORATION, MAYNARD. MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


THIS SOFTWARE IS cyenienee UNDER A LICENSE AND MAY BE USED AND COPIED 
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® 
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® 
® 
$ OT SOFTWARE 0 NY OTHER * 
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® 
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FACILITY: COBOL ARITHMETIC 


“ABSTRACT: 
This module contains the routine which divides two 
quadwords, producing a quadword result. 


VERSION: 1 
HISTORY: 


AUTHOR: 
John Sauter, 26-DEC-78 


MODIFIED BY: 
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-SBTTL HISTORY ; Detailed Current Edit History 


; Edit History for Version 1 of COBDIVG 


- Original from COBOL-74, MODULE FASB4 eOyT ine FAST _DIV_4 

- Make the entry point symbol global. 483 95 AN-1979 ~ 

- Minor editing ¢ oom S. S igts N-1 

> Ego lete rewrite. 519 emor-19 

digit temps and other minor «Aes MLJ 13-Mar-1979 

= Correct DIVP round towarrds zero problem. PDG 12-Jul-1979 

- Make it work corcec eas, wipe overlapping input and output 
operands. R 7 

Cosmetic changes. RKR 21-0CT-79 

Complete rewrite. PDG 15-Jun-1981 

This includes several of the ideas from FAST 

- Added EDIT field for use in checkin's audit eyr*: LB 28-JUL-81 

= Updated copyright date. LB 30-JUL-81 

- Some bug fix. PDG 9-Aug-1981 
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ft «SBTTL DECLARATIONS 
f + INCLUDE FILES: 
75; 
$ 
+ EXTERNAL SYMBOLS: 
0 8 3 NONE 
00 0; 
00 1 
3H mono, 
0 4; NONE 
00 5 3 
000 : 
000 3 
000 8 : PSECT DECLARATIONS: 
00000000 89 .PSECT _COBSCODE PIC, SHR, LONG, EXE, NOWRT 
9000 90 
000 91 ; 
0000 3 : EQUATED SYMBOLS: 
0000 93: NONE 
0000 («94:: 
9000 95 
000 9%: 
0000 97 ; OWN STORAGE: 
0000 (98: 
0000 (99: NONE 
0000 100: 
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-SBTTL COBSDIVQ_R8 


"FUNCTIONAL DESCRIPTION: 
Divides two quadwords, producing a quadword result. 
CALLING SEQUENCE: 
JSB COBSDIVQ_R8 (divisor.rq.r, dividend.rq.r, quotient.waq.r) 
Arguments are passed in R6, R7 and R8. 
INPUT PARAMETERS: 


DIVISOR.rq.r The divi 
DIVIDEND. rq.r The dividend. 


IMPLICIT INPUTS: 
ALL of the trap bits in the PSL are assumed off. 
OUTPUT PARAMETERS: 


IMPLICIT OUTPUTS: 
NONE 
COMPLETION CODES: 
NONE 
SIDE EFFECTS: 
Destroys RO through R8. 
If the divisor equals zero, a ‘‘divide 
divi 
NOTES: 
In comments below, the following conventions are used: 


ed by -1. The result is the largest negative integer. 


Equet ions describe the current state of the registers. 


A primed value represents the an value of a register. 


COBSDIVQ_R8:: 
CMPV $3) #1, (R6), 4(R6) ; Divisor in longword range? 
BNEQ -«- 200 : No, do slow code 


a ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee eB ee ee 
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QUOTIENT.wa.r The result of the division; DIVIDEND/DIVISOR. 


by zero’’ exception is raised. 
antoger overflow occurs if only if the largest negative integer is 


A prefix of ‘'S'' indicates the longword is to be considered signed. 
A prefix of "U"' indicates the longword is to be considered unsigned. 
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; The divisor is a longword, so (hopefully), a single EDIV will work 


Try DIV 

at Ay i forges 

Store quotient 

Store high-order longword 
Return 


(RGR. (R7), RT, RO 


! (RB)+ 
EXTV #31, #1, R1, (RB) 


100$: 


:+ 


oO 


MOVQ (R7), R2 
MOVL (R6), R6 


: We want to compute: 
SR3x2°32 + UR2 = SR6 x (SR1x2°32 + URO) + remainder 


where the remainder is of the same sign as the numerator, 
and is less (in absolute value) than the divisor. 


Grab dividend 
Grab divisor 


o 


ty #31, ahd. RS + Sign extend dividend to R4-R3-R2 
EDIV R6, 7 » R3 :; First EDIV (can't overflow) 


SR4x2°32 + URS" = SR3" = SR6 x SR1 + SR3 
Now try the second EDIV 


§ 
4 
5 
} 
: The divisor is a longword, but a single EDIV doesn't work 
g 
4 
5 
; 
8 
S 
0 
1 
g 
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89 

78 90 DIV R6, R2, RO, RS 

1¢ 31 vc 107$ 
38 We know that the quotient will fit onee 0 3g, bits, but becouse, (in 
94 the true quad result of the divide), 2 doesn't equal bit 31, 
95 EDIV gets overflow and trashes the heoul ts. We borrow (carry?) from 
38 the high Longyore of the quotient so that bits 32 and 31 of the true 
9 qyequore quotient will be Peete — signs of the remainder and the 
38 ividend may differ, but that's okay (for now). 

CD 00 XORLS R6, R3, RS ; Determine sign of true result 

19 1 ss 108$ ; Branch if true cgeylt is negative 

ce $ SUBL2 Re: ee) 3; Subtract SR6x2*°32 from dividend 

0 Z suet Rl ; Increase quotient by 2°52 

co 5 103$: ADDL2 R6, R3 z: Add SR6x2°32 to Cty ideng 

p $ ECL ; Decrease quotient b 2 

B " H's EDIV R6, R2, RO, RS 3; Re-try the second EDIV 

3% 03 : ; SR3x2°32 + UR2 = SR6 x SRO + SRS 
4 11 : SR3*x2*32 + UR2 = SR6_x sri ai 232 + 5R3x2*32 + UR2 
4 i 5 SR6 x SRI x £,3°,”  * 
" 1? : = SR6 x (SR1x2*32 + skO) (getting close) 
18 004 15 cea 108 ; Make RO unsigned 
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51 D7 48 1 DECL R1 
4A 13 108$: 5 
2 18 ; SR3'x2*32 + UR2 = SR6 x (SR1x2*32 + URO) + SRS 
1" 9 ; Fix up the result if the remainder has the wrong sign 
53 54 0=C<‘ OSCE 4A é KORL3 RS, R4, RB 3 Signs same? 
14 1 4 BGEQ 1028 : If so, OK. 
55 D 5 4 TSTL R5 : See if remainder is zero 
10 86 5 5 BEQL 102$ : If so, the quotient is correct 
a3 32 § cstt a i; Was original dividend plus or minus? 
th D $3 8 INCL RO 3; Increment the result 
51 00 D Ben ADWC * R1 
9° 1 5D 0 BRB 1028 
0 D7 OOSF 1 105$: DECL RO ; Decrement the result 
51 9 08 061 ¢ SBwWC #0, R1 
68 0 D 0064 102$: mova RO, (RB) 3 Store result 
05 0067 4 RSB ; And return (remainder = psuedo R5) 
B38 5 200$: 
06 § 3+ 
006 ; 
one HS ; The divisor doesn't fit into a signed longword 
006 40 ;- 
3 4 006 41 CLRL (SP) ; Clear negate flag 
moe D 006A $s MOVQ (R7), R2 ; Grab dividend into R3-R2 
0B 18 006D 4 BGEQ 201$ 3; Skip if positive 
$f 96 O06F 44 INCB (SP) ; Toggle negate flag 
2¢ 5 C 071 45 MNEGL R2, R 3 Negate dividend 
5 99 D 074 46 ADWC #0, R 
53 — 0077 47 MNEGL R3, R 
56 66 «67D «(007A «66248 2018: mova (Rb), RO : Grab divisor into R7-R6 
0B «18 Baoe 49 BGEQ 2028 : Skip if positive 
6— 96 OO07F 50 INCB (SP) ; Toggle negate flag 
56 56 CE 0081 51 MNEGL R6, R ; Negate divisor 
57 29 i) ite 26 ADWC #0, R 
«5 CE 0087 5 MNEGL R7, R7 
008A 54 2028: 
OBA 2? 3; Normalize the divisor 
50 57 01 C9 QOBA 25 BISL3 #1, R7, RO ; Make sure high longword not zéro 
50 50 4E OO8E 58 CVTLF RO, RO : To find most sign} icant bit 
303) 59 ssscece BGTR § ? : Branch if every nine okay 
c 91 60 BLSS 3 : Branch if divisor still negative! 
4 04 0093 61 207$: CLRL R4 needed? ; Clear R4 in case we don't shift 
50 50 08 07 cf 95 6¢ EXTZV #7, #8, RO, RO ; Move the exponent into low byte 
50° 9F 8F 50 83 009A 26 SUBBS RO, #159, RO ; Un-normalize the exponent 
11019 oF 64 BLSS 218 ; needed? : Don't shift right! loses accuracy) 
56 56 50 79 OA} 65 ASHQ = RO, RO, R : Shift the divisor into R7-R6 
O43 $ 3; Shift dividend by same amount, into R4-R3-R2 
21 20 Bh Ag $3 SuBL3 RO, #32, R1 ; Get offset in bits for R4 
54 3 3 1 ai Ag 0 EXTZV. RI, RR. RB. RG 3; Grab portion for R4 (must be extZv!) 
_. = 0 ag 4} 2138 ASH RO, R2, R ; Shift lower two 
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ie 3 Divide R4-R3-R2 by R7-R6 giving R5-R4 remainder R3-R2 
i ; First, divide R4-R3 by R7 giving R4 remainder R3 
5 4 7 ELRL 2 $ Clear high longword Ra. tient 
so 51 53 8? «9B bons 78 ght BF. as, Rt. RO i RESRS / A? = RECRG rem AS 
9 79 sssrecs BVC bs $ na Me if no ates Ray 
He 1D 9 0 ; Branch if we had overflow 
53 0 7D : 1 222$ MOVva RO, R3 3 Quotient in R5-R4, peneinder in R3 
; $ F SR7 x SR4 + SR3 = SR4'x2*32 + URZ' 
See 5 ; Multiply back and subtract to compute remainder by using the formula: 
OBE 3 : SR4'x2°64 + UR3S'x2*32 + Av = (SR7x2°32 + UR6) x SR4 
it 38 : = Rao 3 + UR2 = UR6 x 
50 00 54 56 7A 008 90 295$: émuL RG R4, #0, RO z R1I-RO = UR6 x SRO 
03 26 i EI o¢ 31 BBC A ai’ 269$ > Consider R6 unsigned 
2g 30 2 Been 38 209$:  SUBL RO, R2 : Subtract the product 
5 1 D OCD 94 SBWC R1, R 
gp 18 0000 295 BGEQ §=©=_-20s ; Was estimate too big? | 
4 07 00Dd2 38 204$: DECL RG 3: Yes, decrement estimate 
55 00 D9 00D4 29 SBwC_ 8 #0, RS | 
26 23 D8 Bopt 38 ApOLe RG. : ; And add back divisor 
‘31 Opp 00 ea BLSS 2048 ; And continue until positive 
st 44 8 "  $ See if we must negate the quotient | 
09 8E E9 OODF 04 BLBC (SP)+, 206$ : Skip if no negate required 
2¢ 54 ce Boe B2 ee BB Rd 3; Negate result | 
55 39 CE QOE8 09 MNEGL RS, R5 | 
68 54 7D OOEB 08 206$: mova R4, (RB) ; Store quotient 
05 OOEE 09 RSB ; Return 
ae | 
Ree 1 : This code was moved out of line so that ‘usual’ cases don't cause branches. | 
bore 8 ies, ' 
Boer 1 "> The Stvtoor ~ 2 agile, negative! Tote is because the original divisor 
OEF 1 3 was equal to The goly possible quotients are zero and one, 
+H 18 ; depending on whether the dividend also equals 2°32. 
ee FF OEF 0 éxtzv #31, #1, R3, RS : Get low longuord of result 
68 55 7D OOF4 1 MOoVa = R5,, (RBS ; Store quotient (R6 = 0, remember?) 
8E B2 F7 § TSTL (SP)+ ; Throw away the negate flag 
5 F9 RS 3; And return 
OFA 328 5° 
oA ; : This code was moved out of Line so that ‘usual’ cases don't cause branches. 
OFA += 
OFA § 2918: ; 


=~ 


< 


K 9 
ge Quadwords 15-SEP-1984 $382 :08 penlyes Macro v04-00 Page 8 
66-SEP-1984 10:44:4 COBRTL.SRCIJCOBDIVQ.MAR; 1 (4) 


pet 1vQ_R8 


aon 
22 
RO 
or 
cs 

ze. 


2 SS SS SS SS SS IOS OOOSOSCSOSS <o 


Pee te ee ee ey | = | |) = be in bs en be es Ps Bn ns Pn a Pm a a a ee 


; Since we had overflow, we know that, before normalizing: 


a >= dividend >= 2° 
3; 2°32+7 >= divisor >= 2° 
3 < >= quotient >= 2°31 (true quotient, that is) 
; Also, on overflow, EDIV sets the remaiider to zero. 

+ We can compute the desired results via: 
t € (RG=R3) = R7x2°31 ) / R7 + 2431 = RI rem RO 


; Unless this also couse’ an overflow, in which case, we know that 
; quotient must equal 2°32. 


50 00 57 180000000 of 7A EMUL #14831, R7, #0, RO ; Multiply by -2°31 
0 3 ©6¢O ADDL2 3, RO 3; Compute the difference 
1 0 »d8 ADWC #0, R1 
2) 4 0 ADDL2 R4, RI 
50 «(51 0 7 B IV R7, RO, R1, RO 3; Try the EDIV one more time 
05 1 BVS 292 ; Oh, damn! 
AG 51 #1F OC BBCS #31, R1, 2228 3; Add 2°31 to the quotient 
0 HALT 3; We should NEVER get here 
a. fia. Ae 292$:  SUBL3 R7, R4, RO : Compute the remainder 
51 D4 CLRL R1 3; Set low lLongword of quotient 
55 06 INCL R5 ; Set high longword of quotient 
99 «11 BRB 2228 3; Join common code 
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Symbol table 6- 
COBSD1VO_R8 00000000 RG 01 


tone wen anocmonen nas + 


! Psect synopsis ! 


gn mewn san ane nen en } 


PSECT name Allocation PSECT No. Attributes 
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! Performance indicators ! 


een 2 


Phase Page faults CPU Time Elapsed Time 
Initialization 41 0:00:00.05 88:88:38 85 
ommand processing 138 0:00:00.39 0:00:0 f 
Pass 7 B88: 3: 5 + Se th 
Symbol table sort :00:00.00 H's .00 
Pass 7 BB: 0.46 00:00:02.04 
Symbol table output 1 0:00:00.01 00:00:00.01 
Psect synopsis output 0:00:00.01 0:00:00.01 
Cross-reference output 0:00:00.00 0:00:00.00 
Assembler run totals 32 00:00:01.47 0:00:13.60 


Ihe working set Limit was 1050 pages. 
0 bytes (11 pages) of virtual memory were used to buffer the intermediate code. 
There were 10 pages of symbol table space allocated to hold 1 non-local and 21 local symbols. 
358 source Lines were read in Pass 1, produc ng 8 object records in Pass 2. 

pages of virtual memory were used to define 0 macros. 


Macro Library name Macros defined 
~$255$DUA28: CSYSLIBISTARLET.MLB;2 0 
O GETS were required to define 0 macros. 


There were no errors, warnings or information messages. 


RD NOWRT NOVE 
NOVE 


M* CRO/ENABLE=SUPPRESS10N/D1 SABLE=(GLOBAL , TRACEBACK) /LIS=L1S$:COBDIVQ/0BJ=0BJ$:COBDIVQ MSRC$:COBDIVQ/UPDATE=(ENH$: COBDIVQ) 
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